Limits on spin-independent couplings of WIMP dark matter with a p-type 

point-contact germanium detector 
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We report new limits on spin-independent WIMP-nucleon interaction cross-section using 39.5 kg- 
days of data taken with a p-type point-contact germanium detector of 840 g fiducial mass at the 
Kuo-Sheng Reactor Neutrino Laboratory. Crucial to this study is the understanding of the selection 
procedures and, in particular, the bulk-surface events differentiation at the sub-keV range. The 
signal-retaining and background-rejecting efficiencies were measured with calibration gamma sources 
and a novel n-type point-contact germanium detector. Part of the parameter space in cross-section 
versus WIMP-mass implied by various experiments is probed and excluded. 



PACS numbers: 95.35.+d, 29.40.-n, 98.70.Vc 

About one quarter of the energy density of the universe 
can be attributed to Cold Dark Matter(CDM)[T], whose 
nature and properties are unknown. Weakly Interacting 
Massive Particles (WIMP, denoted by x) are the leading 
candidates for CDM. There are intense efforts to look 
for WIMPs through direct detection of nuclear recoils in 
terrestrial experiments via X-N— elastic scattering. 

Most experimental programs are optimized to study 
WIMPs with mass m x ~10— 100 GeV, motivated by pop- 
ular supersymmetric models. Germanium detectors with 
sub-keV sensitivities were identified and demonstrated 
as possible means to probe the "low-mass" WIMPs with 
m x <10 GeV[2]. This inspired development of p-type 
point-contact germanium detectors (PPCGe) which pro- 
vide such sensitivities with modular mass of kg-scale[3]- 

Our earlier measurements at the Kuo-Sheng Reactor 
Neutrino Laboratory (KSNL, with a shallow depth of 
about 30 meter-water-equivalence) using a 4-element ar- 
ray with a total mass of 20 g and analysis threshold 
of 220 eVee ("ee" denoting electron-equivalence energy 
throughout) have placed constraints on m x >3 GeV[4]. 
The CoGeNT experiment reported data with a 440 g 
PPCGe[S], showing an excess of events at the sub-keV 
range over the background models. A consistent an- 
nual modulation signature was observed. Allowed re- 
gion in the spin-independent j^iV couplings (cr^) was 
derived. Intense interest and theoretical speculations 
in the low-mass WIMP region were generated, some of 



which explored the consequences of additional corrections 
on residual sub-keV surface background^- The low en- 
ergy data of the CDMS and XENON experiments [7j have 
subsequently excluded the allowed region with different 
detector techniques, though these interpretations were 
challenged [8], and the issue remains controversial. It is 
crucial to have independent experiments which can probe 
the CoGeNT allowed region and provide further under- 
standing on the detector response and the nature of the 
sub-keV events in PPCGe. 

We report new results at the sub-keV region with a 
PPCGe having a fiducial mass of 840 g (actual crystal 
mass 926 g) at KSNL. The low-background facilities as 
well as the hardware, trigger and data acquisition (DAQ) 
configurations were described in our previous work [H [5]. 
The PPCGe detector was enclosed by an Nal(Tl) anti- 
Compton (AC) detector and passive shieldings of OFHC 
(oxygen free high conductivity) copper inside a plastic 
bag purged by nitrogen gas evaporated from the liquid 
nitrogen dewar. This set-up was further shielded by, from 
inside out, 5 cm of OFHC copper, 25 cm of boron-loaded 
polyethylene, 5 cm of steel and 15 cm of lead. This 
structure was surrounded by cosmic-ray (CR) veto panels 
made of plastic scintillators read out by photomultipliers. 
Both AC and CR signals are crucial, serving both as ve- 
tos to reject background and as tags to identify samples 
for efficiency measurements. 

Signals from the PPCGe point-contact are supplied 
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FIG. 1: Scatter plot of the PPCGe rise time (logio[r]) versus 
energy. The ro-line corresponds to the BS cut in this analysis, 
with the range for cut-stability test (r-scan) indicated. Typ- 
ical B and S pulses at ~700 eVee are depicted in the insets. 



through a reset preamplifier. The output is distributed 
to a timing amplifier (TA) which keeps the fast rise-time 
information, and to shaping amplifiers at both 6 fxs (SAe) 
and 12 fj,s (SA 12 ) shaping time which provide energy in- 
formation and the triggers. Signals from the PPCGe 
outer surface-electrode are processed with a resistive 
feedback preamplifier and followed by a shaping ampli- 
fier at 4 jus shaping time (SA out ). The TA, SA6,i2,outj 
and AC-Nal(Tl) output were digitized by flash analog- 
to-digital converters at 200 MHz, 60 MHz and 20 MHz, 
respectively. The discriminator and timing outputs of 
the CR panels were also recorded. The trigger thresh- 
old was adjusted such that the trigger efficiency is 100% 
above 300 eVee, verified by test pulser (TST) events. A 
total of 53.8 days of data were taken, where the DAQ 
dead time was 12.6%, measured by random trigger (RT) 
events. Energy calibration was achieved by the internal 
X-ray peaks with the zero-energy defined by the pedestals 
provided by the RT events, while the range in between 
was cross-checked with TST events. 

A cut-based analysis was adopted. There are three 
categories of selection criteria: (i) the "physics ver- 
sus noise events" (PN) cuts differentiate physics sig- 
nals from spurious electronic noise; (ii) the CR and AC 
cuts identify events with activities only at the PPCGe; 
and (iii) the "bulk versus surface events" (BS) cut se- 
lects events at the interior of PPCGe. In addition, 
the efficiencies and suppression factors (e^,Ax) for ev- 
ery selection (A=PN,AC,CV,BS) are measured. They 
correspond to the probabilities of (signal, background) 
events being correctly identified. Only events above 
the electronics noise-edge of 500 eVee were adopted 
in this analysis. The physics samples selected by the 
PN cuts are categorized by "ACV(T)+CRV(T)+B(S)" , 
where ACV(T) and CRV(T) are, respectively, the 



AC- and CR-vetoed(tagged) events which are in anti- 
coincidence(coincidence) with the PPCGe signals, while 
B(S) denote the bulk(surface) samples. The xN candi- 
dates would therefore manifest as ACV+CRV+B events. 

Background suppression with the PN, AC and CR cuts 
and the evaluations of their respective (ex, Ax) follow the 
well-studied procedures as in earlier experiments [41 1^1 1 1 Oj . 
The PN cuts are based on pulse shape characteristics 
and correlations in (TA,SA6,SAi2,SA out ). They suppress 
spurious triggers induced by microphonics effects or the 
tails of pedestal fluctuations. Background induced by 
the preamplifier reset are identified by the timing corre- 
lations with the reset instant. The in situ doubly-tagged 
ACT+CRT events serve as the physics reference sam- 
ples, with which the signal efficiencies epN shown in Fig- 
ure l2c are accurately measured. The majority of the 
electronics-induced events above noise-edge are identi- 
fied (Apjv~1)- The efficiencies for AC and CR selections 
are measured by RT events to be, respectively, eAC>0-99 
and ecR=0.93. The suppressions are Aac=1-0 above the 
Nal(Tl) threshold of 20 keVee, while A C r = 0.92, mea- 
sured by reference cosmic samples in which the energy 
depositions at Nal(Tl) are above 20 MeVee. 

The BS selection, on the other hand, is a unique feature 
to PPCGe. The surface-electrode of PPCGe is a lithium- 
diffused n + layer of mm-scale thickness. Partial charge 
collection in the surface layer (S) gives rise to reduced 
measureable energy and slower rise-time, as compared to 
those in the bulk region (B)[5l [HJ [12]. The thickness 
of the S layer was derived to be /s=1.16 mm, via the 
comparison of simulated and observed intensity ratios of 
7-peaks from a 133 Ba source [T3"]. This gives rise to a 
fiducial mass of 840 g, or a data size of 39.5 kg-days. 

The data are differentiated in B/S event-by-event 
through the rise-time (r) of their TA output. The logio[r] 
versus measured energy ( T) scatter plot and examples 
of typical B/S events are depicted in Figure [T] The 
observed and actual rates are denoted by (B',S') and 
(B,S), respectively. The B&S bands are well separated at 
T>2 keVee, such that (e BS ,A B s)~(100%,100%), or equiv- 
alently, B=B' and S=S'. At lower energy, (B',S') and 
(B,S) are related by the coupled equations: 



B' = e BS ■ B + (1 - A BS ) • S 
S' = (1 - £bs) B + X BS -S 



(1) 



with an additional unitarity constrain: B+S=B'+S'. 

To evaluate (cbs^Abs) at sub-keV energy, it is necessary 
to devise calibration schemes which provide at least two 
measurements of (B',S') where (B,S) are independently 
known. The TST events are inappropriate since their 
TA-output exhibit different pulse shapes from those of 
physics events. Instead three complementary data sam- 
ples, as displayed in Figure^, were adopted: 

(I) Surface-rich events with 7-ray sources — Calibra- 
tions with both high and low energy 7-sources ( 241 Am 
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FIG. 2: The derivation of (e B s,A B s)- (a) Total, B' and B 
spectra from the surface-rich 7-ray ( 241 Am, 137 Cs) and bulk- 
rich cosmic-ray induced neutron calibration data. The rela- 
tive normalizations between (B,B') are fixed by 6bs = 100% at 
T>2.5 keVee. (b) Allowed bands of (£bs,Abs) at threshold 
and at a high energy band, (c) The measured (ebs^bs) and 
epN as a function of energy. 



at 60 keVee and 137 Cs at 662 keVee, respectively) were 
performed to characterize the PPCGe response. As dis- 
played in Figure [2^, the measured B '-spectra are com- 
pared to the reference B derived from full simulation with 
Is =1.16 mm as input. The simulated B-spectra due to 
external 7-sources over a large range of energy are flat 
for T<5 keVee. 

(II) Bulk-rich events with cosmic-ray induced fast neu- 
trons — A 523 g first-of-its-kind n-type point-contact ger- 
manium (NPCGe) detector was constructed. The com- 
ponents and dimensions are identical to those of PPCGe. 
The surface of NPCGe is a p + boron implanted elec- 
trode of sub-micron thickness. There are no anoma- 
lous surface effects in NPCGe. Data were taken un- 
der identical shielding configurations as the PPCGe at 
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FIG. 3: Measured energy spectra of PPCGe, showing the 
raw spectra and those after ACV+CRV cut, as well as +B'. 
The large inset shows the (eBS,A B s)-corrected ACV+CRV+B 
spectrum, with a flat background and L-shell X-ray peaks 
overlaid. The small inset depicts the residual spectrum su- 
perimposed with that of an allowed (excluded) cross-section 
at m x =7 GeV. 

KSNL. The ACV+CRT condition selects cosmic-ray in- 
duced fast neutron events without associated 7-activities, 
which manifest mostly(~85%) as bulk events in PPCGe. 
Accordingly, as depicted in Figure [2^,, the measured 
ACV+CRT data with NPCGe is taken as the B-reference 
and compared with those of ACV+CRT+B' in PPCGe. 

Using calibration data (I) and (II), (cbs^bs) are evalu- 
ated by solving the coupled equations in Eq. [2] As exam- 
ples, the three allowed bands on (ebSjAbs) at threshold 
and at a high energy band are illustrated in Figure [2|d. 
The different orientations of the bands are consequences 
of the different depth distributions of the samples, which 
give rise to different B:S ratios. The bands have common 
overlap regions, indicating the results are insensitive to 
the event locations. The surface-rich 7-events and the 
bulk-rich cosmic-ray induced neutron-events play com- 
plementary roles in constraining Abs and 6bs, respec- 
tively. The results are depicted in Figure [2};, with epN 
overlaid. Consistency is checked by comparing the mea- 
sured strength of the in situ Ga-L X-ray peak at 1.3 keVee 
after BS-selection to that predicted by the corresponding 
K-peak at 10.37 keVee. 

The raw spectra as well as those of ACV+CRV and 
ACV+CRV+B' are depicted in Figure [3j The peaks 
correspond to known K-shell X-rays emissions from 
the cosmogenically-activated isotopes. The (cbSjAbs)- 
corrected spectrum of ACV+CRV+B is shown in the 
large inset. Errors above ~800 eV are dominated by sta- 
tistical uncertainties, while those below are defined by the 
systematics from the calibration procedures. When the 
expected L-peaks and a flat 7-ray background are sub- 
tracted, the residual spectrum corresponds to %N candi- 
date events and is shown in the small inset. In addition, 
the stability of the results is studied over the variations 
of tq within the r-scan range of Figure [I] The corre- 
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TABLE I: The p-values for the stability hypothesis on 
the key parameters following a r-scan of Figure [I] The 
ACV+CRV+B rates are insensitive of the choice of to. 



Energy Range (keVee) 


0.5-0.7 


1.5-1.9 


ess 


< 10 _b 


< 10~ 5 




< 10" 5 


< 10~ 5 


ACV+CRV+B' 


0.24 


0.16 


ACV+CRV+S' 


0.17 
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ACV+CRV+B 


0.57 
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FIG. 4: Exclusion plot of spin-independent x~N coupling, 
superimposed with the results from other benchmark experi- 
ments. 

sponding changes in (e S s,As,gjB',S',B) are derived. The 
p- values for the stability hypothesis are listed in Table [I] 
The ACV+CRV+B rates are stable and insensitive to the 
choice of To, even with (ess, Ass) exhibiting significant 
variations in the expected directions. This demonstrates 
that the calibration procedures are valid and robust. 

Constraints on are derived via standard approach 
with conventional astrophysical modelspQ (local density 
of 0.3 GeV/cc and Maxwellian velocity distribution with 
v =220 km/s and u esc =544 km/s). The rates of WIMP- 
induced events due to spin-independent interaction can- 
not be larger than the residual spectrum. The quench- 
ing function in Ge is evaluated with the TRIM software 
which matches well with existing dataJTU]. As illustra- 
tion, xN recoil spectrum due to an allowed (excluded) 
cr^ at m x =7 GeV is shown in the inset of Figure p| 
Exclusion plot of a^ N versus m x is displayed in Fig- 
ure |4j Bounds from other benchmark experiments are 
superimposed [51 [71 H"4"]. An order of magnitude improve- 
ment is achieved compared to our previous results [3]. 
Part of the DAMA, CRESST-II and CoGeNT allowed 
regions are probed and excluded. We note that an excess 
remains in the sub-keV region not yet accounted for in 
this analysis, the understanding of which is the theme of 
our on-going investigations. 



Studies continue on PPCGe and NPCGe at KSNL. 
Projects on improvement of electronics and sub-noise- 
edge analysis [12] are being pursued. Dedicated dark mat- 
ter experiment CDEX with sub-keV germanium detec- 
tors are taking data at the new China Jinping Under- 
ground Laboratory (C JPL) [T3] . This facility provides the 
attractive features such as a rock overburden exceeding 
2400 m and horizontal drive-in access. 
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